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We present a relativistic calculation of the saturation properties of nuclear matter which contains 
the correlation energy. Pion loops are incorporated on top of a relativistic Hartree-Fock (RHF) 
approach based on a chiral theory. It includes the effect of nucleon structure through its response 
to the background chiral invariant scalar field. All the parameters which enter the RHF calculation 
are fixed or strongly constrained by hadron phenomenology or lattice data. The new input for the 
correlation energy is the Landau-Migdal parameter g' governing the short-range part of the spin- 
isospin interaction. We find that the inclusion of the correlation energy improves the description of 
the saturation properties of nuclear matter. 
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I. INTRODUCTION 

In a set of recent papers [l|, [3, we have proposed a relativistic description of nuclear matter where the nuclear 
binding is obtained with a chiral invariant background field governing the medium modifications of the non-pionic 
piece of the chiral condensate. As a starting point wc take the point of view that the effective theory has to be 
formulated in terms of the fields associated with the fluctuations of the chiral quark condensate parametrized in a 
matrix form, W = a + if ■ tt. The dynamics of these fields is governed by an effective potential V{a, t?) having a 
typical mcxican hat shape reflecting a broken (chiral) symmetry of the QCD vacuum. Such an effective theory can 
be seen as emerging from a low (space-like) momentum expansion of a bozonized Nambu-Jona-Lasinio action and the 
connection of the scalar " sigma" field a with physical scalar mesons is not necessarily implied [3| • Explicit application 
to the description of nuclear matter of such a N JL model has been performed in ref . Q . 

As proposed in a previous paper Q an alternative and very convenient formulation of the resulting sigma model is 
obtained by going from cartesian to polar coordinates i.e., going from a linear to a non linear representation, according 
to: W = a + if ■ n = SU = (/^r + s) exp (ir • 0-,^/ f-n)- The new pion field corresponds to an orthoradial soft 
mode which is automatically massless (in the absence of explicit chiral symmetry breaking) since it is associated with 
rotations on the chiral circle without cost of energy. The new sigma meson field S*, which is a chiral invariant, describes 
a radial mode associated with the fluctuations of the "chiral radius" around its vacuum expectation value /tt. This 
expectation value plays the role of a chiral order parameter around the minimum of the effective potential and the 
medium can be seen as a shifted vacuum. With increasing density, its fluctuation s — S — f.^ governs the progressive 
shrinking of the chiral circle and the evolution of the nucleon mass. Our main physical assumption proposed in ref. Q 
consists in identifying this chiral invariant s field with the sigma meson of nuclear physics and relativistic theories of 
the Walecka type, or, said differently, with the background attractive scalar field at the origin of the nuclear binding. 
This also gives a plausible answer to the long-standing problem of the chiral status of Walecka theories 0] • 

It is nevertheless well known that in such chiral theories, independently of the details of the modelling, tadpole 
diagrams associated with the mexican hat potential automatically generate attractive three-body forces destroying 
saturation Q . The origin of this failure can be attributed to the neglect of the effect of nucleon substructure linked 
to the confinement mechanism as already pointed out in some of our previous works [U, 0, • It was demonstrated that 
nuclear matter stability can be restored once the scalar response of the nucleon depending on the quark confinement 
mechanism is properly incorporated in a way inspired from the QMC model [lo| . In particular it was shown that 
a sufficiently large scalar nucleon response supported by lattice data [T]| is able to generate three-body repulsive 
forces overcompensating the attractive ones coming from chiral tadpoles @. This theoretical framework has been 
implemented in nuclear matter calculation at the Hartee level [l| . In a subsequent work we also incorporated non 
relativistically the pion loop correlation energy. A full relativistic Hartree-Fock (RHF) calculation was then done in 
Q allowing to reproduce also the asymmetry properties of nuclear matter. The aim of this paper is to combine these 
results and to provide a fully relativistic calculation of the pion loop correlation energy on top of the RHF calculation. 

Our article is organized as follows. The second section is a brief summary of our previous works and gives the 
main results of the RHF approach. Section 3 is devoted to the formalism that we utilized for the calculation of 
the polarization propagators entering the correlation energy. Finally in section 4 we make the connection of the 
present work with our previous non relativistic calculation of the pion loop energy and numerical results are given 
and discussed. 
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II. SUMMARY OF THE RHF APPROACH 



111 our previous paper we used the following lagrangian which includes the effect of the scalar field s, omega, rho 
and pion exchanges 



C = ^-27^9^* + Cs + + Cp + (1) 



with 



£s = - V{s) + ^df'sd^s 



The form of Mn{s) 



Mn{s) - Mn + .95 s + i KNs + ^ ) (3) 



reflects the internal nucleon structure through the presence of the scalar response of the nucleon, kats, which can 
be constrained by lattice data analysis According to this analysis we obtained an estimate of the dimensionless 
parameter: C = (/^/2Mjv) k^vs — 1.25. This parameter governs the contribution to the in-medium nucleon mass 
and generates three-body repulsive forces which are necessary to get nuclear matter stability and saturation. The 
scalar mass parameter mo- = 800 MeV is also taken from the lattice data analysis. In addition the nucleon mass may 
very well have higher order derivatives with respect to the scalar field. In practice, as in our previous works [J, we 
also introduce a cubic term. Hence the scalar susceptibility becomes density dependent : 

HNS{S) = „ y = KNS 1 + "T W 



and vanishes at full restoration, s — — where s is the expectation value of the s field. Hidden in the above 
Lagrangian is the explicit chiral symmetry breaking piece 

^xSB =ca = - ^Tr{U + s) cxp {if ■ (p„ / f^) ~ cs - ^ (5) 

which generates the pion mass term with the identification c = fT^rn^. It is thus implicit that neglecting the higher 
order terms in the exponent, the self-interactions of the pions are omitted. Notice that the only meson having a 
self-interacting potential V{s) is the scalar meson s. We take it in practice as in the linear sigma model with the 
inclusion of the explicit chiral symmetry breaking piece : 
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" 2 " + 2/. ' + 8/2 • 

The other parameters {g^^, gp, Kp, qa and the known meson masses) will be fixed as most as possible by hadron 
phenomenology. In particular we use the Vector Dominance picture (VDM) which implies the identification of Kp 
with the anomalous part of the isovector magnetic moment of the nucleon, i.e., Kp = 3.7. The omega meson should 
also possess a tensor coupling but, according to VDM the corresponding anomalous isoscalar magnetic moment is 
~ 0.13. Since it is very small we neglect it here. 
The energy density of symmetric nuclear matter writes : 

^ — ^kin+Hartree ~t- ^Fock (7) 

with 

= y (k . - + M.i-s) - j + Virs) + - j (8) 



3 



The equation of motion for the classical scalar field s is : 

_ v^s + V'is) = ~g*s PS with g*s = ^^^^ = gs + knss + .■■ (9) 

The effective nucleon mass M* and the effective momentum k* come from the Hartree-Fock equations. As discussed 
in rcf. @] in the non relativistic limit the orbital wave functions arc imposed; these are simply non rclativistic plane 
waves. Consequently if the system is not too relativistic, the results of the calculations should not depend very 
much on the choice of the wave functions. Said differently, the results obtained with another basis than the fully 
self-consistent Hartree-Fock (HF) basis would deviate from the HF results only by tiny relativistic effects as shown 
in 3 . For this reason we introduce the Hartree basis which ignores the Fock term in the nucleon self-energy. In that 
case the effective mass M* coincides with Mn{s), the effective momentum k* coincides with the normal momentum 
k and the effective nucleon energy becomes E*{s) = y^p^ + . One advantage will be the strong simplification 

of the calculation of the polarization propagator entering the calculation of the pion loop correlation energy. 

The Fock contribution from scalar exchange is given by 

^Fock -—J (27r)3 (27r)3 (k - k')^ + mf 

where to*^ ~ V"{s) + kns Ps is the effective scalar mass governing the propagation of the scalar fluctuating field. 
The other Fock term contibutions (omega, rho and pion) to the energy density are listed in the appendix B of ref. 0] . 
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III. FORMALISM FOR THE CORRELATION ENERGY 



For the calculation of the pion loop correlation energy we generalize to the relativistic case the approach used in 
0. The correlation energy associated with pion loops is calculated using the well known charging formula method: 

i^r^^V^e;oop = 3V^ f^'-^ J (^/x^'^;M.(^'q;^) (nr(^,q;A) - n(;:(.;,q;A)) . (11) 

Va;pj/(ci;, q; A) is the residual interaction in the axial- vector channel wich contains together with the pion exchange 
potential the effect of Migdal short-range correlations (Landau-Migdal parameter g') taken in a covariant form {q^ = 
Lu^ -cf) according to [ll[ll]: 

In our calculation we will choose a dipole form factor v{q) with a cutoff A = 0.98 GeV, such that the resulting 

contribution to the free nucleon sigma term, cr^j^\ is cr^'' = 21.5 MeV, in agreement with previous determinations 
[U, [l^, [l^ . In the charging formula all the coupling constants at the interaction vertices are systematically weighted 
with a A factor. This is the origin of the A^ appearing explicitly in the residual interaction. n^''(aj, q; A) is the 
full polarization propagator in the axial-vector channel in presence of the interactions weighted by A. It is defined 
according to : 

nr(g) ^yj dit^t')dvdr'e^'^<^-^'^ (0|(-^)T(^(^7V^) (x), (^tV^) (^')) |0) ■ (12) 

Notice that we are actually dealing with correlators in an isovector channel. We omit for simplicity the isospin Pauli 
matrices Tj since the only effect in symmetric nuclear matter is to multiply at the very end the result for the correlation 
energy by a factor three. To get the genuine correlation energy, the first order term, i.e., the mean field one, has to 
be substracted since it is already incorporated as the Fock pion exchange term in presence of short range correlations 
Q . The explicit expression for these mean- field fl" polarization propagators will be given below. 

In the non relativistic limit this correlation energy involves polarization propagators (or interactions) only in the 
spin-isospin channel: the pion exchange interaction is of pure longitudinal spin-isospin nature whereas the short 
range interaction contains both longitudinal and transverse pieces 01 ■ However, non relativistically, the transverse 
spin-isospin polarization propagator is also affected by rho meson exchange. Hence to keep this connection with the 
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non relativistic case 0], we also incorporate the rho meson exchange which is already present in our previous RHF 
calculation. The rho exchange contributes to the correlation energy according to: 

E';'"p^ve,,o,^3v f^-^ J (0i/T^p^M.('^'q;^) (nr(^'q;^) - n(;;(c.,q;A)) . (13) 

Vp^^{uj, q; A) is the rho meson exchange interaction: 

q; A) = - X' gl ^ v^q). 

For the vector coupling constant we take as in our previous paper the VDM value: gp = 2.65. The rho meson 
propagator should also contain a q^q'^ term but it can be omitted since current conservation implies qfjlp'^ = 0. 
Notice that we take for simplicity the same form factor, i.e., the same momentum cutoff, as in the pion exhange. We 
will discuss this point at the end of the paper and show that taking a harder form factor does not significantly alter 
the conclusions. The vector-isovector polarization propagator is defined (omitting again the isospin Pauli matrices) 
according to: 

nr(9) = dit-t')drdr'e^^<-~^'^ (0| T (^r^(q)V') (x), {^r^''iq)^P) [x')^ |0) . (14) 
The pNN vertices are given by: 

For the tensor coupling we again take the VDM value, Kp = 3.7. As previously discussed the axial and vector meson 
correlators mix (in the non relativistic limit they both contain a transverse spin- isospin piece). Furthermore we also 
have to incorporate the mixed polarization propagators : 

n;^:(<z) = y d(t-t')drdr'e^ {0\ {-t)T (^{i:T^^{q)^) {x), (^7VV^) (a^')) |0) (15) 

Kpiq) = ^ j d{t-t')dvdv'e^^<^-^"^ (0|(-z)t(^(^7V^) (x), (V;rt-(q)^) (x')) |0). (16) 

The polarization propagators are calculated using a RPA scheme. They are solution of the coupled Dyson equations : 

nr('z) = k:{<i) + K:{i)yaM<i)K''{i) + K:p{i)Vp:o.MKa^<i) 

= Klaii) + n(;,"Jg)K;„;3(g)nf''(g) + n'^';{q)Vp.MiHa{q) 
nr(9) = n(;;(g) + n':,;[q)Vp,^p[q)Kii) + no7a('?)K;„/3(g)nf;(g) 

where the A dependence in the arguments have been omitted for simplicity. The mean-field polarization propagators 
entering the above Dyson equations are: 

= ^.tvD [G(p)7VG(p-t-g)7V] (18) 

KU^) = -2 / ^,trn [ G{p) T^^ G{p + q) Yl"] (19) 

K:p{<l) = -2 / i0^,i^D [ G{p) Gip + q) Ft''] (20) 

n^;(g) = -2 1 ^trn [ G{p) F^^ G{p + q) Vl^] . (21) 
The mean-field nucleon propagator has the form : 

p'^ — M*'^ + irj 
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where A'p is the occupation number and it is understood that the pure vacuum piece of the mean-field polarization 
propagators is dropped. 

To solve the RPA problem and to get a compact form for the correlation energy, we will use a projector technics 
p^ . [l7| . For this purpose we introduce the following four- vectors and tensors : 

,7^ = (1,0), = = ^23) 

q- fj^ 

The tensors L, R and T satisfy projector properties and are mutually orthognal. In addition they also satisfy the 
normalization conditions, i^'^L^^ = R^'^ R^i, = 1 and T'^T^^ = 2. The residual interaction can be decomposed on 
these covariant tensors: 

/ \ 2 r 

{q^D^iq) - g') i"^ - g' T^^ - g' R^^^ 



= VaL + VaT T"" + VaR R^" (25) 

yril) = -^■^ gl^'^M -r^ — o {L"'' + T^"" + R"") = Vp (L'"' + T'"' + i?'^") . 

The polarization propagators entering the expression of the correlation energy are actually the various projections on 
the above tensors, namely: 



iiaLiq) = L^, (g), TiMq) = ttT^. nr(<z), n.^Xg) = r^, (g), 



UpLiq)^L^,Il^/iq), UpTiq)^^T^uU^p''iq), Upii{q) = R^^U'^'^ (q). (26) 
The mean-field polarization bubbles have similar projections which means that they can be decomposed according to: 

U^iq) = L^^'-UoaLiq) + T^^'-noaLiq) + R^'-noaB.iq) 

= ^''^ nopL(g) + R^'' nopfl (g). (27) 

We see on these expressions that the rho channel polarization propagator has no projection on the longitudinal channel 
which remains true for the full propagator. There is an additionnal mean- field bubble which mixes the axial and rho 
channels. It has the following tensorial structure: 

KUq) = K:p{i) = Tr no6(g) with = ^f'""'fip<i^ 

T6^.r'''^ = 2, Tg^^Tea.^n, TrT:=T^\ (28) 

The antisymmetric Tg tensor is actually orthogonal to the L and R tensors. Consequently the mixing beween the 
axial and rho channel proceeds only through the transverse channel associated with the tensor T. We also introduce 
the full axial-rho polarization propagators projected on the Tg tensor 

^pa{q) = ^pa■,^.u{q) ^ap{q) = ^Tg^"^ Tlap,p,M). (29) 
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The correlation energy can be expressed in terms of the projected polarization propagators: 



j^loop j^loop 



E. 



loop 



e: 



loop 



E 



loop 

aR ~ 



E 



loop 



pR 



aL 

sv 

3V 
SV 
3V 



+°° idm 
-oo 27r 



'OO 
+ 00 



— OQ 



00 



i dcv 
2^ 



j-iloop 

dq 

(2^ 

dq 
dq 

dq 

(27r)3 

dq 

[2-, 



TPIOO} 

EpR 
1 d\ 
T 



KL(^,q;A) (nQL(tj,q;A) - noQL(w,q;A)) 



/■I 

—2VaT{^,q;X) {UaT{iO,q;X) - UoaT{io,q; X)) 
—2VpT{oJ,q;X) {IlpT{oJ,q;X) - nopT(^,q;A)) 

/■I 

J — Ki?(w,q;A) {UaR{uj,q;X) - UoaR{uj , q; X)) 



r^ldlO f dq f'dX u m ^ TT ^ 

J ~2^ J (2^ J (npfl(a',q;A) - nop_R(tj, q; A)) 



(30) 



The Dyson equations for the projected polarization propagators can be obtained from the original ones (eq. I17[) by 
projecting them on the various channels L, T or R. After straightforward algebraic manipulations one obtains: 



n„L 


= HoaL 


+ IIoqL VaL HaL 




= HoaJ?, 


+ Hoai?. VaR Ilafl 




= HoaT 


+ HoaT VaT ^aT + Hoe l^p ^paT 




= Hoe + 


Hoe V'aT HaT + HopT V^p HpaT 


n^T 




+ Hp Vp HpT + Hoe Kt ^apT 




= + 


Hoe Vp IIpT + HoaT 14t HapT 


^pR 




+ nop_R Vp TlaR - 



(31) 



From the solution we get the particular combinations entering the expression of the correlation energy: 

VaL IIoaL 



VaL HaL 
VaR ^aR 
VaT^aT + VpWpT 
VpUpR 



1 — VaL IIoaL 

VaR. HoaiJ 
1 — VaR^OaR 

KtHq^t + VpHopT +2KTVp(n^6 - np^TnopT) 

(1 - VaT^OaT){l - VpU^pr) " Kt l^p Hgg 
Vpliop 



1 - K,n, 



p ^^OpR 



(32) 



The integration over the varying coupling constant A is immediate and we get the final expression for the various 
contributions to the correlation energy: 



jTjloop _ j^. 



e; 



loop 
aL 



j^loop 
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loop 
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dq 
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ln(l - VailloaL) + VailloaL 



i duj 
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i duo 
~2^ 



dq 

(2^ 



In [(1 - UoprVp) (1 - UoaTVaT) - VaTTioeVpUoe] + UoprVp + HoaTKr 

ln(l - VaRUoaR) + VaRHoaR (uJ , q) 

ln(l - VprHopr) + VprHqpr (w,q). 



(^^,q) 



(33) 



The leading term is obtained by expanding the log to second order in VHq which yields the result of second order 
perturbation theory which is always negative according to a basic result of quantum mechanics. In particular the 
second order pion loop is embedded in the axial-longitudinal contribution to the correlation energy It can be directly 
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compared with the iterated pion exchange (the so-called planar diagram) appearing in medium chiral perturbation 
theory. If the ChiPT calculation is regularized with a cutoff [l^, a negative result is also obtained but much larger 
than in our approach where the pion exchange is strongly screened by short-range correlations. 

The analytical structure of the polarization propagators considered as a fonction of w is such that it has a cut on 
the real axis and these propagators are analytic in the first and third quadrants. In other words the continuous set 
of poles lies below the real axis for lo positive and above the real axis for lo negative, the same as for the pion and 
rho propagators. For this reason the (practical) calculation can be done using a Wick rotation. Each four momentum 
integral will be calculated according to: 

In practice as an input of the calculation all what we need is the real part of the bare polarization propagators to make 
an analytical continuation to uj = iz . The explicit form of these polarization propagators are given in the appendix. 



IV. DISCUSSION AND NUMERICAL RESULTS 

a. First order polarization propagators and non relativistic limit. The various mean-field projected propagators 
can be written generically as: 

no,(g) 



with 

KoaLiq;p) 
Koariq; p) 
KoaRiq; p) 

KopTiq;p) 
KapRiq; p) 

where k* — Kp(iV/*/A/jv) appears as an in- medium modified rho tensor coupling constant and = p — p • qq, with 

p|, = p^ — (p • q)^, is the transverse component of the nucleon momentum. 

The difference between relativistic and non relativistic calculations has essentially three distinct origins: the first 
is the incorporation of nucleon-antinucleon excitations, the second originates from the use of Dirac spinors which 
translates into the structure of the Kj functions and the third is pure kinematics (for instance Ep replaced by 
M* + p^ /2A1*). To appreciate the role of the NN excitations we rewrite the denominator appearing in eq. [35las: 

KM 1 ^ K,{q) / 1 1 

Ep + 2p ■ q + if] 2 Ep Ep+q - Ep+q + Ep + iij lu + i?p+q + Ep - ir] 

The first term represents the forward going p — h bubble and the second term represents the backward going NN 
bubble. For the energy denominator associated with q changed in — g, we write similarly: 

Kj{-q) 1 ^ _ K,{-q) / 1 1 

Ep q^ - 2p ■ q + irj 2EpEp-q \UJ + i?p-q - Ep - irj lj - -Bp-q - Ep + irj 

The first term represents the backward going p — h bubble and the second term represents the forward going NN 
bubble. The inclusion of antinucleons is at best questionable since there are many excitations far below the NN 
threshold which should be incorporated before. In addition the internal consistency would require to take into 



4dp fNj 



I 



K,{q;p) 



(27r)3 \Ep q'^ + 2p-q + ii] 



+ (q^ -q) 



^dp o._^P^P+q 



(27r) 



Er, 



Kj{q;p)S{q^ + 2p-q)eiEp + lo) 



(35) 



2M*^ 

2 M*^ + p-q + Pt 
-2p^ 

I + K* 



2M* 



2M*{\q\ Ep - ujp-q) 



l±^)\»2M-.p.,4 + pJl-(^)"/ 



= -2M* 



2AI* 



q^ {p-q + 2p|) - 2p| + kU^. 



(36) 
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account vacuum polarization effects but the restructuring of the QCD vacuum has probably little to do with virtual 
NN excitations. However at the formal level they are essential to maintain the covariance of the calculation and the 
simple tensorial structure and symmetry properties (g^HJ^" = 0). For that reason we decided to keep them in our 
calculation as many authors [13, [3- Moreover their practical effect is a small 0{1/AIn) relativistic correction. 

We now come to the various Kj functions (eq. I36[) . The axial longitudinal bubble IIoql (pion channel) is actually 
combined with VaL which contains the pion exchange whith an explicit —q^ = — cj^ factor. This is a potentially im- 
portant difference with the non relativistic approach where the pion coupling contains only (pure p-wave coupling). 
However since the dominant energy range is jIM* (quasi-elastic peak) the relativistic correction induced at 

the level of the vertex is of order ^^(l/M^). The conventionnal non relativistic limit is obtained by ignoring the 
antinucleon terms and by replacing l/2i5p i5p_|_q by l/2Af*^ so that Kj{q)/EpE-pj^q ~ 1 with the result : 

4dp _ „ f 1 



<aL = -no(c^,q)^- 



(2^)^ 



p+q 



En 



irj 



E, 



p+q 



En 



IT] 



(39) 



We also see that the K functions for the transverse axial channel and the rho-R channel arc governed by the dominant 
non relativistic term (±2M*^), the terms in and being of 0(l/il/^) order. The term involving p ■ q = EpOj — p q 
can also be seen of second order since uj ~ q^/2M*. The analysis of the rho transverse channel is a little particular 
since the first term is proportionnal to q^ and apparently of order 0{1/M^). However the factor 1 -I- k* is about 5. 
Hence for typical momentum |q| ~ 300 MeV the factor (1 -I- k*)/2 Af* |q| turns out to be of order unity such that 
the transverse rho and axial-rho bubbles survive in the non relativistic limit. We thus conclude this discussion by 
summarizing the bare bubbles surviving in the non relativistic limit: 
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(40) 



In the non relativistic limit, the equations (j32p reduce to: 
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1 - V^i^^Ho 



1 



1 - VpUo 

where the non relativistic longitudinal and transverse spin-isospin interactions are given by: 



(41) 
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Consequently, in this limit the correlation energy writes: 
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(43) 



The first two terms represent the longitudinal and transverse spin-isospin contributions as in ref. 0. The last 
term has nothing to do with spin-isospin physics but corresponds to the contribution of the time component of the 
vector interaction which is anyway very small {—0.65 MeV at normal nuclear matter density). As a byproduct we 
can estimate the contribution to the correlation energy of sigma and omega exchanges that we have not considered 
explicitly in this paper. This omission can be justified with the following argument. In the non relativistic limit there 
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FIG. 1: Binding energy of nuclear matter with g' = 0.7. The values of the parameters and the coordinates of the saturation 
point are given in Table IT] The full line corresponds to the full result, the dotted line to the mean-field (RHF) contribution 
and the dashed line to the correlation energy. 

is no difference between the scalar density and the vector density. The relevant interaction appearing in the RPA 
summation for the a + uj channel is: 



Since the sigma, rho and omega have very similar masses around 800 MeV and since the leading term to the correlation 
energy involve two meson exchanges lines, one can estimate 



where the 1/3 is a trivial isospin factor. Taking = 8.1, gg = 6, one obtains i?g°f^ ~ —0.8 MeV which is also very 
small and hence justifies the use of the mean-field approximation for the omega and sigma exchanges. 

b. Numerical results and discussion. As explained before the spirit of the approach is to study nuclear matter 
properties with various parameters fixed as much as possible by hadron phenomenology and lattice data. We only 
allow a fine tuning for the ujNN coupling constant, g^^, around the VDM/quark model value, 3gp ~ 7.95 and for the 
nucleon scalar response parameter around the value estimated from the lattice, Ciatt — 1.25. The only nuclear physics 
input is the Landau-Migdal parameter for which we first take the vastly used value g' = 0.7 compatible with the 
most recent data analysis [20j . For g^ = 7.6 and C = 1.14 wc obtain and excellent reproduction of nuclear saturation 
properties as shown on fig. [TJ As compared with our previous pure Hartree-Fock work Q the C parameter is reduced. 
Another satisfactory result is the value of the compressibility modulus K = 9p^d^e/dp^ = 250 MeV, smaller than our 
previous values and very close to the accepted value around 240 MeV. This can be understood by looking at figure [1] 
where the mean-field and the correlation energy contributions to the binding energy per nucleon are separately shown. 
The important point is that the correlation energy displays a non linear behaviour with density which helps saturation 
to occur: a smaller scalar response parameter is needed hence reducing the incompressibility. It is also interesting 
to examinate the various pieces contributing to the correlation energy. We see on the last column of table [Til that, 
with the chosen value of g' , there is an almost complete screening of the pion exchange and the correlation energy 
is dominated by the transverse channel. At the saturation point we find for the longitunal axial piece E^°^^ = —0.6 
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FIG. 2: Same as figure 1 but for g' = 0.5. The values of the parameters and the coordinates of the saturation point are given 
in Table U 



g' 


0.7 


0.5 




7.6 


7.3 


G 


1.14 


1.3 


P/PO 


1.00 


1.00 


E/A (MeV) 
K (MeV) 


-15.97 
250 


15.87 
270 



TABLE I: Values of the parameters and coordinates of the saturation point when the correlation energy is included on top of 
the RHF calculation. 

MeV whereas the dominant transverse piece is = -5.5 MeV. We also see that the axial R channel which is a 

pure relativistic correction is negUg ible: e'-°°p = -QMMeV. Finally the relativistic correction to the rho exchange 
energy which is not of spin-isospin nature is E''"^^ = —0.65 MeV. 

It is also interesting to discuss the influence of some inputs of the calculation although they arc constrained by 
accepted phenomenology. If the g' parameter is decreased to the value g' = 0.5 the saturation properties can be 
also reproduced (sec fig El) but the compressibilty modulus is larger: K ~ 270 MeV (see Table H]). In that case the 
screening of the pion exchange is less pronounced E^°^^ ~ —2.4 MeV but the transverse piece is reduced e}^°p = -2.1 
MeV. Actually this transverse contribution exhibits a minimum around g' = 0.2 which corresponds to the strongest 
compensation between attractive rho exchange and the contact interaction. This can be seen on Tabic El where the 
various contributions to the correlation energie are given at the normal nuclear matter density for various values of 
g' . An important point is that decreasing g' brings nuclear matter close to pion condensation [2l|. The onset of pion 
condensation corresponds to a pole at zero energy in the pion propagator or in the axial longitudinal polarization 
propagator at some critical momentum q^, namely (1 — VaL^i:iaL){'^ ~ ^iQc) = 0. This is illustrated on the left 
panel of fig. [3] where E^^^^ is displayed for various values of g'. For the case g' = 0.15 the calculation ends at a 
density p = 1.8 po which is the critical density for pion condensation. It turns out that the relativistic calculation 
disfavors pion condensation as compared with the non relativistic one. For instance in the non relativistic calculation 
pion condensation occurs at p = 2.2 po for g' = 0.3 whereas it is beyond 3 po in the relativistic case. Close to 
pion condensation the longitudinal correlation energy strongly increases which is reminiscent of a critical opalescence 
phenomena [2^ . 



-20k 



0.5 



1.5 
P/P, 



2.5 



-20k 



0.5 







1.5 

P/Po 



2.5 



FIG. 3: Longitudinal axial contribution to the correlation energy (left panel) and total correlation energy (right panel) as a 
function of density for various values of g' . Solid curves: g' = 0.7; dashed curves: g' = 0.5: dotted curves: g' = 0.3; dot-dashed 
curves: g' = 0.15. 





g' = o 


g' = 0.1 


g' = 0.2 


g' = 0.3 


g' = 0.4 


g' = 0.5 


g' = 0.6 


g' = 0.7 


rploop 


X 


12.5 


-8.8 


-6.0 


-4.0 


-2.4 


-1.3 


-0.6 


^loop 


-2.6 


-1.5 


-0.9 


-0.9 


-1.4 


-2.1 


-3.7 


-5.5 


rploop 
^aL 


-0.00 


-0.00 


-0.00 


-0.01 


-0.01 


-0.02 


-0.03 


-0.04 


rploop 


-0.65 


-0.65 


-0.65 


-0.65 


-0.65 


-0.65 


-0.65 


-0.65 



TABLE IL The various contributions (MeV) to the correlation energy at normal nuclear matter density for various values of 
g' with C = 1.14 and g^i = 7.6. The cross for the longitudinal contribution for the g' = case means that pion condensation 
occurs. 

The magnitude of the correlation energy is controlled by the cutoffs entering the form factors. For the t:N form 
factor, it is fixed by the pion cloud contribution to the pion-nucleon sigma term. However there is no real reason 
to keep for the rho exchange contribution the same form factor. Therefore we have changed the cutoff entering the 
pN form factor from = 980 MeV to Ap = 1500 MeV. It turns out that the transverse contribution is not much 
increased and the saturation curve can be obtained with only a tiny change of the C and g^^ parameters, the main 
effect being the increase by a factor three of the relativistic correction, E^°l^ to the rho meson exchange. This is 
illustrated in Table [1111 

However even if this form factor may seem to be quite soft we also have to keep in mind that a more careful treatment 
of short range correlation generates a momentum dependence of the g' parameter together with a tensor h' parameter 
in such a way that the residual interactions in both channels vanish at high momentum (Beg- Agassi-Gal theorem) 
hence accelarating the convergence of the momentum integration p3| . Another point is the neglect of delta- hole 
bubbles which contribute significantly to the correlation energy according to our previous non relativistic calculation 
. However we found in 0] that it gives a contribution almost exactly linear in density which is phenomenologically 
undistinguishable from omega exchange. In practice a very small increase of g^o is expected to simulate this delta-hole 
bubble contribution. 



V. CONCLUSION 



In this work we have developped a fully relativistic RPA framework for the calculation of the correlation energy. 
We have included the effect of pion loops, as well as the rho ones, on top of a RHF calculation. Our description 
of saturation properties is almost parameter free. The nuclear binding is ensured by a background chiral invariant 
scalar field associated with the radial fluctuations of the quark condensate. In order to reach saturation we have 
incorporated the response of the embedded nuclcon to the background scalar field. It generates three-body repulsive 
forces which allows saturation to occur. All the parameters entering the Hartree-Fock description, including the scalar 
nucleon response, have been fixed or constrained by hadron or QCD phenomenology. 
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Ap [MeV] 


980 


1500 


1500 




7.6 


7.6 


7.8 


C 


1.14 


1.14 


1.15 


K [MeV] 


250 




270 


E'""" [MeV] 


-6.8 


-9.4 


-9.4 


Kr [MeV] 


-0.6 


-1.6 


-1.6 



TABLE III; Dependence of the correlation energy with the cutofT entering the pN form factor. The first column corresponds 
to the original calculation with g' = 0.7 and the second column gives the modification of the correlation energy when changing 
the cutoff. The last column shows the new parameters when the saturation point is readjusted. 

The Landau-Migdal parameter g' entering the correlation energy is taken from the phenomenology of spin-isospin 
physics. We have found that the magnitude of this correlation energy is moderate, of the order —10 MeV per nucleon, 
due to the strong screening of the pion exchange by short range correlations. Eventhough this contribution is not very 
large its non linear behaviour with density improves the result of the RHF mean-field calculation for the description 
of bulk properties of nuclear matter. This confirms our previous conclusion obtained in a non relativistic evaluation 
of the correlation effects 0] ■ This result seems to be robust in the sense that this conclusion is not altered if we vary 
the pN form factor or the g' parameter around g' = 0.7, the standard value. 
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APPENDIX: Explicit form of the bare polarization propagators after 
Wick rotation 



In this appendix we give the expressions for the various mean-field polarization propagators once the integration 
over the angle between the transferred moment q and the nucleon momentum p has been performed. They are given 
after the Wick rotation (eq. [34|) which is used to perform the practical calculation of the correlation energy. We use 
the following tool functions inspired from [2^ : 

"(92-2|p||q)2+4zX2\ 



lni(z, |q|) = In 

2:ln2(2;, |q|) = —2zatan 
The expressions of the various propagators are : 



(<z2 + 2|p||q|)2 + 4z2i?2y 
&\p\h\zEp 



(g2)2 _4p2q2 +4^2^2 



HoaT 
IIoaL 



1 f pdp 

1 1 

1 1 

1 q2 
47r2 |q|3 



+ lap 

|p||q| ^+4i?pzln2(z,|q|; 



/ ^ {(-4Ar2^2 ^ 4|p|2^2 ^ (^2)2) ^ 8|p||q|^2 ^ ^Ep^x^^{z, |q|)} 

|^4Ar2lni(z,|q|) 



pdp 

En 



^El - z2 + [q|2 + 4i^M*2 ) lni(z, |q|) + 8|p||q|. 



-z^ + |qp 



-4£;pzln2(2;, |q|) 



<|qplni(z,|q|)+ ^ ^ 



{(-4iir2|q|2 + 4i?2g2 ^ (^2)2) ^ 8|p||q|g2 + 4£;pzln2(z, |q|)} 



n 



OpR 



1 (j2 [pdp 
47r2 |q|' 



2 V2A/jv, 

{ {^El + q^) lni(z, |q|)8|p||q| + AEpZ l^z, |q|)} - 2K;\qf lni(z, |q|) 



-* „2 



{ (z2|q|2 + 4p2q2 + isy) lni(z, |q|) + 8|p||q|z2 - Aq^EpZ In^iz, |q|)} 



06 



I ^{2i?,lni(z,lq|) + .ln2(z,lq|)} 



(44) 



[1] G. Chanfray and M. Ericson, EPJA 25 (2005) 151. 

[2] G. Chanfray and M. Ericson, Phys. Rev C75 (2007) 015206. 

[3] E. Massot and G. Chanfray, Phys. Rev. C78 (2008) 015204. 

[4] L.S. Celenza, Huangsheng Wang and CM. Shakin, Phys. Rev. C763 (2001) 025209. 

[5] W. Bentz and A.W. Thomas, Nucl. Phys. A696 (2001) 138. 

[6] G. Chanfray, M. Ericson and P.A.M. Guichon, Phys. Rev C63 (2001)055202. 

[7] B.D. Serot, J.D. Walecka, Adv. Nucl. Phys. 16 (1986) 1; Int. J. Mod. Phys. E16 (1997) 15. 

[8] A.K. Kerman and L.D. Miller in "Second High Energy Heavy Ion Summer Study", LBL-3675, 1974. 

[9] M. Ericson and G. Chanfray, EPJA 34 (2007) 215. 

[10] P.A.M. Guichon, Phys. Lett. B200 (1988) 235. 

[11] D.B. Leinweber, A.W. Thomas and R.D. Young, Phys. Rev. Lett 92 (2004) 242002. 

[12] M. Nakano et al, Int. J. Mod. Phys. ElO (200) 459. 

[13] M.F.M. Lutz, Phys. Lett. B552 (2003) 159. 

[14] I. Jameson, A.W. Thomas and G. Chanfray, J. Phys. G18, (1992) L159 . 

[15] M.C. Birse and J.E. McGovern, Phys. Lett. B292 (1992) 242. 

[16] A. W. Thomas, P. A. M. Guichon, D. B. Leinweber and R. D. Young Progr. Theor. Phys. Suppl. 156, 124 (2004); 
|nucl-th/041101 4| 

[17] L.SrCelenza, A. Pantziris and CM. Shakin, Phys. Rev. C45 (19992) 205. 

[18] N. Kaiser, S. Fritsch and W. Weise, Nucl. Phys. A697 (2002) 82. 



14 



[19] J.C. Caillon and J. Labarsouque, Phys. Rev. C59 (1999) 1090. 

[20] M. Ichimura, H. Sakai and T. Wakasa, Progr. Part. Nucl. Phys. 56 (2006) 446. 

[21] A.B. Migdal, Rev. Mod. Phys. 50 (1978) 107. 

[22] J. Delorme and M. Ericson, Phys. Lett. B76 (1978) 182. 

[23] G. Chanfray, J. Delorme and M. Ericson, P hys. Rev . C31 (1985) 1582. 

[24] L.P. Leinson and A. Perez, Ar.XiV nuc l-th/0307025| 



